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The output pulses of a commercial high-power femtosecond fiber laser or amplifier are typically 
around 300-500 fs with a wavelength around 1030 nm and 10s of /iJ pulse energy. Here we present 
a numerical study of cascaded quadratic soliton compression of such pulses in LiNbC>3 using a type 
I phase matching configuration. We find that because of competing cubic material nonlinearities 
compression can only occur in the nonstationary regime, where group- velocity mismatch induced 
Raman-like nonlocal effects prevent compression to below 100 fs. However, the strong group velocity 
dispersion implies that the pulses can achieve moderate compression to sub- 130 fs duration in 
available crystal lengths. Most of the pulse energy is conserved because the compression is moderate. 
The effects of diffraction and spatial walk-off is addressed, and in particular the latter could become 
an issue when compressing in such long crystals (around 10 cm long) . We finally show that the second 
harmonic contains a short pulse locked to the pump and a long multi-ps red-shifted detrimental 
component. The latter is caused by the nonlocal effects in the nonstationary regime, but because it 
is strongly red-shifted to a position that can be predicted, we show that it can be removed using a 
bandpass filter, leaving a sub-100 fs visible component at A = 515 nm with excellent pulse quality. 

PACS numbers: 42.65.Re, 42.65.Ky, 05.45.Yv, 42.70.Mp, 42.65.Hw, 42.65.Jx, 42.65.Jx 



I. INTRODUCTION 

Pulsed fiber laser systems are currently undergoing a 
rapid development, and by employing the chirped pulse 
amplification (CPA) technique high-energy femtosecond 
pulses can be generated with /iJ-sub-mJ pulse energies 
1] . Combined with the fact that the fiber laser technology 
offers a rugged, cheap and compact platform, ultrafast 
fiber CPA (fCPA) systems could compete with solid-state 
amplifier systems. However, the gain bandwidth of the 
Yb-doped fibers typically used for lasing in the 1.0 /zm 
region is considerably lower than competing solid-state 
materials (such as TkSapphire crystals). Thus, due to 
the build up of an excessive nonlinear phase shift Yb- 
bascd fCPA lasers are often limited to a pulse duration 
that typically is sub-ps at best (around 500 — 700 fs) for 
~ 100 fiJ pulses 2] while shorter pulses can be reached 
(~ 250 fs) for - 30 //J pulses @. 

Efficient external compression methods are therefore 
needed. A prototypical compressor consists of a piece of 
nonlinear material, where a broadening of the pulse band- 
width occurs by self-phase modulation (SPM), followed 
by a dispersive element (gratings or chirped mirrors) that 
provides temporal compression. With this method (using 
a short piece of fiber as nonlinear material) 27 fs sub-//J 
pulses were generated from 270 fs 0.8 /j.3 pulses from an 
fCPA system Alternative methods consist of using 
long (0.5 m or more) gas cells or filaments [f| as non- 
linear material, and this works with pulse energies from 



50 /.tJ to around 1 mJ (limited in part by self-focusing 
effects) or possibly even higher energies @. 

Using soliton compression both the SPM-induced pulse 
broadening and dispersion-induced compression occur in 
the same material [7[. However, as self- focusing solitons 
require anomalous dispersion this can only be achieved 
in the near-lR through strong waveguide dispersion. This 
means using specially designed fibers, such as micro- 
structured fibers. Fibers have a very limited maximum 
pulse energy of a few nJ, albeit large mode-area micro- 
structured solid-core and hollow-core fiber compressors 
can support up to 1 /iJ 

Unfortunately the pulse energy from fCPA systems 
lies exactly in the gap between these methods. We will 
here study a compression method that can compensate 
for this. It is a soliton compressor based on cascaded 
quadratic nonlinearities [0l— Fill] . see Fig. [1] This has sev- 
eral advantages: As it relics on a sclf-dcfocusing nonlin- 
earity, there are no problems with self-focusing effects, 
and multi-mJ pulse-energies can be compressed. More- 
over, solitons require normal instead of anomalous dis- 
persion, implying that solitons can be generated in the 
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FIG. 1. (Color online) The cascaded quadratic soliton com- 
pressor studied here: the Yb fiber laser produces energetic 
longer pulses 100 fs) that are launched collimated in a 
quadratic nonlinear lithium niobate crystal, where the phase- 
mismatched type I SHG process compresses the input pulse. 
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visible and near-IR. Finally, it is extremely simple as it 
relies on just a small piece of quadratic nonlinear crystal, 
preceded only by a lens or a beam expander [L?! ]. 

The basis for the cascaded quadratic soliton com- 
pressor (CQSC) is phase-mismatched second-harmonic 
generation (SHG). The cascaded energy transfer from 
the pump (fundamental wave, FW) to the second har- 
monic (SH) and back imposes a strong SPM-like non- 
linear phase shift on the FW, whose sign can be made 
self-defocusing [H, [TJ]. Thereby the FW pulse can be 
compressed with normal dispersion , and soliton com- 
pression becomes possible in the visible and near-IR [l(| • 

In this paper we investigate the CQSC in a type I 
lithium niobate (LiNbC>3, LN) crystal, where the goal 
is to perform moderate compression of longer fs pulses 
from fCPA systems at the Yb gain wavelength of 1030 
nm. We show that in order to overcome the detrimen- 
tal cubic nonlinearities the phase mismatch has to be 
chosen so low so that the compression occurs in the so- 
called nonstationary regime. This regime is dominated 
by group-velocity mismatch (GVM) effects, and exactly 
the large GVM is a well-known drawback of using LN 
in the near-IR for SHG. However, when only moderate 
compression is desired, the soliton order can be kept low, 
and we show through numerical simulations that reason- 
able pulse quality can be achieved and that up to 80% 
of the pulse energy is retained in the central spike. The 
compression limit is found to be around 120 fs FWHM, 
which is a limit set by the GVM effects. The compression 
occurs in a crystal of reasonable length, 10 cm. This is 
possible only because LN has a very large 2. order dis- 
persion. Finally, we show that bandpass filtering of the 
SH actually can lead to a very clean sub-100 fs visible 
pulse with around 0.1% conversion efficiency. 

In this paper we first discuss the general compression 
properties of LN in a cascaded type I SHG interaction 
setup in Sec. [TTJ and then show some numerical simu- 
lations in Sec. IIIII of pulses coming from two different 
commercially available fCPA systems. We conclude in 
Sec. IIVI The properties of LN are discussed in App. 
and App. [B] discusses the anisotropic Kerr nonlinear re- 
sponse of LN. Appendix [Cl and iDl discuss the conversion 
relations between Gaussian and SI units for cubic non- 
linear coefficients and Miller's rule, respectively. 

II. TYPE I COMPRESSION PROPERTIES OF 
LITHIUM NIOBATE CRYSTALS 

With the CQSC high-energy few-cycle compressed 
pulses can be generated, as was experimentally observed 
at 1250 nm [151 ]. However, the first studies performed at 
800 nm were plagued by GVM effects, that prevented 
reaching the few-cycle regime || [ill EH ■ These studies 
used a /3-barium-borate (BBO) crystal in a type I SHG 
00 — > e configuration, where the FW (ordinary polar- 
ization) is orthogonal to the SH (extraordinary polariza- 
tion) and where birefringent phase matching is possible 



by angle-tuning the crystal. BBO is in many respects an 
ideal nonlinear crystal: it has low dispersion, a very large 
transparency window, and a reasonably strong quadratic 
nonlinearity relative to the detrimental cubic one. As we 
have shown in previous theoretical and numerical studies, 
BBO provides an excellent compression of longer pulses 
to ultra-short duration at the Yb gain wavelengths f]~6l — 
[l8J. The problem with BBO is that good quality waveg- 
uides are not supported and that it is very difficult to 
grow long crystals. Especially the latter is important if 
only moderate compression of longer pulses is desired. In 
moderate soliton compression most of the pulse energy is 
conserved in the compressed pulse, and the pulse has a 
reduced pedestal. The problem is that compression will 
only occur after a long propagation length. 

We therefore turn here to LN, which is a widely used 
quadratic nonlinear crystal for IR frequency conversion. 
LN is attractive due to extremely large effective quadratic 
nonlinearities (up to 10 times larger than BBO), that can 
be accessed through a quasi-phase matched (QPM) type 

SHG phase matching configuration where FW and SH 
have identical polarization. However, here we study LN 
in a type I configuration as BBO. The effective quadratic 
nonlinearity is more than twice as large as in BBO. 

LN is usually not considered very suitable for SHG 
of short pulses in the near-IR because the SH becomes 
very dispersive; thus, the FW and SH group velocities 
are very different resulting in large GVM. This is also 
why LN has not been used in the near-IR as nonlin- 
ear medium for the CQSC, for which GVM is a very 
detrimental effect. Another disadvantage for the CQSC 
is that the Kerr nonlinear response is several times larger 
than BBO, which counteracts the advantage of the large 
quadratic nonlinearity of LN. Therefore the CQSC exper- 
iments done so far using LN were done in the telecommu- 
nication band and exploited QPM in a type configura- 
tion fl9| , where effective quadratic nonlinearity is around 
three times larger than what can be achieved in a type 

1 configuration. However, we now show that type I LN 
offers a quite decent compression performance without 
having to custom design a QPM grating. 

A. Solitons with cascaded quadratic nonlinearities 

In cascaded quadratic interaction the FW effectively 
experiences a Kerr-like nonlinear refractive index. This 
is in addition to the cubic (Kerr) nonlinearities that are 
always present in all media. We can write the total re- 
fractive index of the FW [see Eq. (|C2|) ] 

n = rix + \ \£\\ 2 n cnh - lc = ni + Iin^ ubic (1) 

where n\ is the FW linear refractive index, £\ is the 
FW electric field, and I\ the FW intensity. It is typi- 
cal to report the nonlinear refractive index relative to 
the electric field, n cu bi c , or to the intensity, n* nhic . We 
have here for simplicity neglected cross-phase modula- 
tion (XPM) contributions since they are small in cas- 
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caded SHG. As mentioned we have contributions from 
both cascaded quadratic and cubic Kerr nonlinearities 



"■cubic — n SHG 



l Kcrr,ll 



(2) 



where ?%crr 11 i s t ne SPM Kerr nonlinear refractive index 
of the FW (see App. [B]for details on the notation etc.). 
The contribution from the cascaded quadratic nonlinear- 
ities can in the large phase mismatch limit (AkL 1, 
where L is the crystal length) be approximated as [ill 



'SHG 



c£oAinfrt2A/c 



(3) 



where d c g is the effective \^ nonlinearity. For Ak = 
k 2 — 2fci > the cascaded contribution is negative, i.e. 
self-defocusing. Here kj = ^tt/Xj is the wavenumber. 

The effective quadratic nonlinearity of the type I oo — > 
e interaction for the 3m crystal class (LN, BBO) is 



d e ff = sin — d^2 cos sin 3<fi 



(4) 



where the angles are defined in Fig.[5]in App.lBl Choosing 
6 = —ir/2 gives maximum nonlinearity (see App. |A")) . 

In cascaded quadratic soliton compression the aim is to 
get ng HG < and |ng HG | > n^ orr u as to achieve a total 
sclf-defocusing cubic nonlinearity. The soliton interaction 
can then be described by an effective soliton order [ItI ] 
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where iV SH G = L D ikxI h 



'SHG 



is the soliton order of 



the self-defocusing cascaded quadratic nonlinearity, and 



N K 



Kerr, 11 



is the soliton order of the 



material Kerr self-focusing cubic nonlinearity. The FW 
dispersion length is Ld,i = ^n/|&± |> where k^ is 
the FW group- velocity dispersion (GVD). We generally 
use the following notation for the dispersion parameters 
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FIG. 2. (Color online) Properties at Ai = 1.03 /im when angle- 
tuning the LN crystal: (a) Phase mismatch, (b) GVM param- 
eter, (c) GVD of FW (red) and SH (blue), and (d) the spatial 
walk-off angle p. The effective quadratic nonlinearity neglect- 
ing (black) and including (dashed red) spatial walk-off are 
shown in (e) and (f) is the total cubic Kerr nonlinearity ((2} 
from cascaded quadratic nonlinearities and Kerr SPM (using 
18 x 1(T 20 m 2 /W, see App.©. 



i 

^Kerr 



anomalous and self-defocusing solitons are no longer sup- 
ported. The SH GVD (blue) is about 3 times larger than 
the FW GVD. 

Since the type I critical phase matching is employed, 
the walk-off angle p — arctan[tan(0)rig/rig] — (valid for 
a negative uniaxial crystal) is nonzero, see Fig. Hfd). In 
Fig- IH e ) h is apparent that d e g is largely unaffected by 
walk-off. However, walk-off does set a limit to the effective 
interaction length between the pump and the SH as we 
will discuss later. 



B. Linear and nonlinear response of LN at 1.03 /im 

Selecting Ai = 1.03 /im, the operating wavelength of 
most Yb-based fiber laser amplifiers, the properties of 
LN are summarized in Fig. [2] the phase mismatch (a) 
becomes small at 9 ~ 1.3 radians (70 — 75°). As shown 
in (e) in this range cZ e ff — 5.2 pm/V, and the total non- 
linear refractive index (f), as expressed by Eq. can 
become negative, implying that the cascaded nonlinear- 
ity is stronger than the Kerr nonlinearity. This happens 
for Afc < 62 mm- 1 (or 6 > 70.4°). At 9 = 75.8° phase 
matching is achieved, after which «g HG > and thus 
sclf-focusing. 

GVM is very large, see Fig.^b), which as we will see 
later sets a strong limitation to the compression perfor- 
mance. The GVD is shown in Fig.[2jc), and importantly 
FW GVD (red) is large and normal (i.e. positive). It will 
stay normal until Ai > 1.9 /im, after which it becomes 



C. Compression diagram for type I LN 

We now generalize to other wavelengths and summa- 
rize the type I compression performance of LN in Fig. [30. 
This compression diagram shows the different compres- 
sion regimes for the CQSC as the wavelength and the 
phase mismatch is varied. 

Above the red curve the total nonlinear refractive index 
is focusing n^ ubic > 0, so solitons are not supported since 



The specific crystal chosen in this work is 1% MgO doped stoi- 
chiometric LN, as the MgO doping gives a much higher material 
damage threshold. Also 5% MgO doped congruent LN would 
work well. See App. [S] for more details about the crystal. 
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FIG. 3. (Color online) Compression diagram for 1% MgOisLN 
at room temperature and aligned for type I SHG. For vari- 
ous pump wavelengths Ai the choice of phase-mismatch pa- 
rameter Ak affects the compression. In order to excite soli- 
tons the phase-mismatch must be kept below the red line 
(Ak < Afc Cjmax ), because otherwise the material cubic nonlin- 
earities are too strong (nKerr,n > I^-shgD- Optimal compres- 
sion occurs when the cascaded nonlinearities dominate over 
GVM effects (Ak > Ak sr , above the black line). We have also 
indicated the operation wavelengths of Yb and Er doped fiber 
lasers. The red line uses Miller's rule to estimate the nonlin- 
ear quadratic and cubic susceptibilities at other wavelengths, 
cf. Eqs. l|DT ]l -l[D2 )l . and uses n' KcllM = 20 x 10" 20 m 2 /W for 
A = 0.78 /im (see App.[B]for an extended discussion). 



The optimal compression occurs in the so-called "com- 
pression window" [16], where the soliton compressor 
works most efficiently because solitons are supported in 
the stationary regime. The diagram shows a compression 
window for type I LN in the regime Ai = 1.6 — 1.9 /im. 
Unfortunately in this range there are no fCPA systems. 

Fortunately, as we will show also in the nonstationary 
regime compression is possible, as long as the effective 
soliton order is low enough. This is what we will try to 
exploit in the regime around Ai ~ 1.03 — 1.06 /im. 

Coming back to Ai = 1.03 /im we observe that solitons 
are supported for when <C Afc < Afc Cjmax = 62 mm -1 . 
However, when getting too close to Afc C:inax the inten- 
sities required to observe solitons become very large im- 
plying excessive Kerr XPM effects and increased Raman- 
like GVM effects [l8j]. On the other hand for Afc too 
small the cascading limit ceases to hold, and also the 
compressor performance decreases due to excessive GVM 
effects [18]. In fact, as a rule of thumb the compres- 
sion limit in the nonstationary regime (in which the 
system will always be for Afc ~ 0) the compression 
limit is roughly given by the pulse duration for which 
^coh = Lgvm 7 where L co h = 7r/|Afc| is the coherence 
length and Lqvm = Af so iiton/Mi2| is the dynamic GVM 
length of a sech-shaped soliton. With "dynamic" we mean 
that the GVM length changes as the soliton compresses. 
Thus, in the nonstationary regime the limit is 
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FWHM 
limit 
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lAfcl 



(8) 



the FW GVD is normal. The curve is found by setting 
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Below the black curve the compression performance is 
dominated by GVM effects (nonstationary regime) while 
above it is dominated by cascaded effects (stationary 
regime) . The curve is to second order given by [l6j 



Afc Rr — 



d 2 

"12 



2k 



(2) 



(7) 



where dyi = k^ — fcj is the GVM parameter and kif' is 
the SH GVD. The lower this curve is the better because 
this implies that the chance of observing solitons in the 
stationary regime increases. Thus, the very large GVM 
parameter di2 is detrimental because it pushes the curve 
upwards. Instead the huge SH GVD values, see Fig. [2] 
(c), are actually helping to push the curve downwards. 
Therefore a large SH GVD can actually be beneficial for 
clean soliton compression. 



,(2) 



2 A more accurate transition can easily be calculated numerically 
using the full SH dispersion operator [T^l , which we have done 
in what follows. 



where the factor in front of the fraction is the conversion 
factor to FWHM for a sech-shaped pulse. Obviously as 
Afc approaches the phase matching point the soliton can- 
not compress to short durations. We numerically found 
the optimal compression point in the Afc = 35 — 50 mm -1 
regime, and with the best results for Afc = 45 mm , for 
which rigjjQ = 25 x 10~ 20 m 2 /W. 



D. Predicting the compression performance 

The next step is to estimate what the compression per- 
formance could look like. Here the scaling laws come 
into the picture, which can be used to predict the prop- 
agation distance for optimal compression z out , the com- 
pression factor f c and the pulse quality Q c 17). 

As we have pointed recently |l6[, it is the phase mis- 
match and the GVM (zero and first order dispersion) 
that really control the compression properties. The only 
requirement to the second order dispersion is that FW 



3 Note that this expression differs with a factor of 7r/2 from the 
limit Tji gjjG = 2|di2/Afc| that we suggested in [T^l ; this is 
purely an empirical choice. 

4 Note that the scaling laws presented here are only ball-park fig- 
ures when used in the nonstationary regime as they were found 
in the stationary regime. 
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(2) 

GVD is normal k\ > as to support solitons. Other- 
wise as we discuss below the FW GVD is basically just 
determining the optimum compression length. The SH 
GVD instead plays a minor role in the compression prop- 
erties, cf. Eq. 0. Our initial idea was to exploit that LN 
is quite dispersive when pumped at Ai ~ 1.0 /im, so the 
very large FW GVD makes it possible to compress the 
pulse in a short crystal. 

So why and when is it interesting to increase GVD as 
to compress in a short crystal? Obviously, the crystals 
have length limits, which for LN is around 100 mm. The 
optimal compression point scales as [l7| 



Zopt 



0.44 2.56 



N 3 „ 

err 



- 0.002. 



(9) 



where zo = §£d,i is the soliton length [20(. So the point 
where the pulse compression is optimal depends on the 
effective soliton order, the input pulse duration and the 
FW GVD. Therefore since quality LN crystals are maxi- 
mum 100 mm long, the CQSC works best when the soli- 
ton order is large and the GVD length is short. But when 
the soliton order is large, the detrimental effects due to 
GVM are strongly increased [H, Ell , m particular in the 
nonstationary regime. Therefore, in the case we study 
here clean compression can only be done with low soli- 
ton order, and therefore the FW GVD must be large as 
to ensure compression in realistic crystal lengths. 

A downside to the large GVD is the following: given 
that some effective soliton order is required then since 



N e g oc IinLu.i oc Ti n y Ii n /\k[ 2 ^\ we have that a large 
GVD gives a short GVD length, and thus larger inten- 
sities are needed to excite a soliton. The same problem 
is found for short input pulses, say from a Ti:Sapphire 
amplifier. However, this is only an issue if operating 
with intensities close to the damage threshold, which is 
not the case here: the intensities are moderate (L m <C 
100 GW/cm 2 ), and instead our issue is to get the soli- 
tons to compress in a crystal that is not too long. 

The compression factor f c = T- ln / Ai opt , where Ai opt 
is the pulse compressed pulse duration at z pt > 

is also 

affected by the effective soliton order [TtJ 



f c = 4.7(V cff - 0.86) 



(10) 



The pulse quality can also be predicted, and is defined 
as the ratio between the compressed pulse fluence with 
that of the input pulse. It scales as [l7| 



Q c = [0.24(V cff - 1) 



l.n 



(11) 



We can use this to calculate the compressed pulse peak 
intensity J opt = Q c fJin and energy E opt = Q c E in . An 
advantage of using low soliton orders is that Q c remains 
high, and thus the compressed pulse retains most of the 
initial pulse energy. 



E. Compression performance of fCPA systems 

Let us use these scaling laws to predict the compression 
performance of fCPA systems. High-energy femtosecond 
pulses from fCPA systems use both Yb doped and Er 
doped gain fibers. Since fCPA systems are diode pumped 
with a wavelength just below 1.0 /im the quantum effi- 
ciency of Yb doped systems is higher, and therefore the 
majority of commercial and scientific systems prefer to 
use Yb over Er. Most systems operate at the A ~ 1.03 /im 
Yb emission line and can for low pulse energies (< 15 /iJ) 
generate pulses as short as 250 fs, while higher pulse ener- 
gies result in longer pulses (currently 50 /iJ 450 fs pulses 
is the state-of-the-art for commercial systems). In Er am- 
plifier systems much lower pulse energies are available, 
typically 1 — 3 /iJ and 500 — 700 fs pulses at A = 1.55 /im; 
such low pulse energies and long pulse duration mean 
that only very low soliton orders can be excited, and thus 
the CQSC can only achieve very moderate compression 
occurring in very long crystals. 

The basis for the following case studies and numerical 
simulations is therefore a couple of commercially avail- 
able Yb-based fCPA systems, both operating at 1030 nm. 
Case (1) is a Clark MXR Impulse giving 15 /iJ 250 fs 
FWHM pulses, which represents a system giving quite 
short, yet still reasonably energetic pulses as a starting 
point. Case (2) is an Amplitude Systemes Tangerine □ 
giving 50 /iJ 450 fs FWHM pulses, which represents a 
system with more energetic but also longer pulses. 

The two cases are studied together taking Afc = 
45 mm -1 . Figure Ufa) shows that in case (1) we need 
to focus the pulses to wo < 600 /im to observe solitons: 
in this regime Fig. BJd) shows that the Rayleigh length 
zji = ttw^/X is only 5-6 times larger than the optimal 
compression point z opt of around 100 mm. This is border- 
line at the risk of experiencing diffraction problems. Even 
increasing or decreasing the waist does not improve this 
ratio much. In case (2) instead, the increased pulse energy 
makes solitons appear already at wq ~ 1.6 mm, despite 
the longer pulse duration. This means that diffraction 
should be less of an issue: in Fig.[4](d) the pulse compres- 
sion point relative to the Rayleigh length of the focused 
beam is significantly smaller in case (2). 

Fig.HKc) indicates that the spatial walk-off in the crys- 
tal can become an issue: the crystal should be shorter 
than the spatial walk-off length L wo = wq/ tan/5 ~ wo/p 
to ensure proper interaction between the FW and the SH, 
but evidently the pulse compression lengths in both cases 
are at least a factor of 2-3 longer than the spatial walk- 
off length. It might therefore be necessary to compensate 
for this by using two crystals, one inverted relative to the 
other so the walk-off direction in the 2. crystal is inverted 
with respect to the 1. crystal [2l| . 



5 http: / / www. dark- m xr~com| 

6 http:/ /www. amplitude-systemes.com 
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FIG. 5. (Color online) Numerical simulation of soliton com- 
pression in LN with Ai = 1.03 fj,m, lf WHM = 250 fs, 
Ait = 45 mm -1 and 7V cft = 1.4 (implying I- m = 6.9 GW/cm 2 ). 
The FW pulse shown in (a) compresses to At opt = 126 fs 
(FWHM) after propagating 91 mm. The SH time plot (c) and 
FW (b) and SH (d) spectra are also shown on a logarithmic 
scale, and Uj are normalized to the peak input FW electric 
field. In (e) and (f) cuts are shown at the optimal compression 
point z = 91 mm (corresponding to the white line in the 2D 
plots). Note that the SH in (e) is magnified 100 times. 



III. NUMERICAL SIMULATIONS 



FIG. 4. (Color online) Practical operation range of the LN 
type I compression system at Ai = 1.03 /im for Afc = 
45 mm -1 . The plots show the predicted behaviour when the 
FW waist wo is varied. The two cases are (1) pump pulses 
with T^ WHM = 250 fs and 15 /xj pulse energy, and (2) pump 
pulses with r£ WHM = 450 fs and 50 /iJ pulse energy. The 
curves in (b)-(f) are calculated based on N e g shown in (a) by 
using the scaling laws [l7| that hold for N e g > 1. 



An alternative solution to the walk-off problem is to 
turn to a noncritical phase matching scheme, where p — 
0. This happens for 9 = or n/2, see Fig.[2ld). Of course 
this removes the possibility of tuning the phase matching 
via 9, and one has to turn to temperature tuning of Afc. 
The temperature needed to get to the desired operation 
point (Afc ~ 40 — 50 mm -1 ) can be estimated using the 
temperature dependent Sellmeier equations [22I ] , and our 
calculations indicate that it should happen already at a 
temperature of around 45° C. This would make an easy 
solution to the walk-off problem. 

The strong GVM implies that compression of Yb-based 
systems can only occur in the nonstationary regime, see 
Fig. [3J Thus, unless N c g is close to unity the GVM in- 
duced Raman-like effects dominate, and the FW pulse be- 
comes extremely distorted and very poorly compressed. 
Actually, as a rule of thumb it never makes sense to use 
N e g larger than what is sufficient to reach the limit ex- 
pressed by Eq. ©, and typically even an N c g smaller 
than that. The limit is drawn as a dotted line in Fig. life), 
and it is reached around wq = 400 /jm in case (1) and 
wo = 800 /xm in case 2. 

Finally, Fig. \4%b) shows that quite moderate input in- 
tensities must be used to achieve solitons in both cases. 
This is related to the quite long input pulse durations. 
Furthermore, Fig. 0Jf) shows that the low soliton orders 
conserve most of the pulse energy in both cases. 



We here present numerical simulations of the two cases 
using a plane-wave temporal model based on the slowly 
evolving wave equation (see more details in [17] and refer- 
ences therein), which includes self-steepening effects and 
higher-order dispersion. This model is justified as long 
as diffraction is minimal, which we assume is the case 
when the crystal length is much shorter than the Rayleigh 
length, and when spatial walk-off is minimal. This re- 
quirement will be discussed further below. 

A. Case (1): 250 fs 15 /iJ pulses 

For the 250 fs 15 /iJ pulses from a Clark laser system 
we found that the best compression was obtained with 
N c s ~ 1.3 — 1.5. This soliton order can be achieved with 
15 [i J pulse energy when the pump is focused to around 
w = 400 (j,m, see Fig.HJa). 

The theoretical compression factor for such soliton or- 
ders is f c = 2-3, i.e., a At opt ~ 80-125 fs FWHM com- 
pressed pulse is predicted. In Fig. [5] we show the results 
of a simulation with N c s = 1.4. This soliton order gave 
the best compression: a slightly asymmetric Ai opt = 126 
fs (FWHM) pulse is observed after 91 mm of propaga- 
tion, see (a) and cut in (e). The compression is not quite 
as strong as predicted by the scaling law (flQ|) . but this is 



7 




J3 0.5 



FIG. 6. (Color online) Simulations as in Fig. [5] but with in- 
creasing N e g. The red curve corresponds to the optimal com- 
pression point from Fig. [5je), while the black curves show 
what happens as the effective soliton order increases (making 
the optimal compression point occurring sooner). 



because the scaling laws are based on pulse compression 
in the stationary regime. On the other hand the pulse 
quality is large, Q c — 0.82, so most of the pulse energy 
is retained in the central compressed part, and the pulse 
pedestal is also very small. These are the main advan- 
tages of soliton compression with low soliton orders. 

The FW spectrum (b) experiences upon propagation 
SPM-like broadening, where the blue-shifted shoulder 
clearly dominates; this is a sign of the cascaded quadratic 
nonlinearities dominating, and the fact that it is blue 
shifted is related to the negative sign of d%2. 

In the SH time-plot (c) we observe the strong GVM 
first inducing a weak component quickly escaping from 
the central part of the pulse, and later the GVM induces 
the characteristic DC-like trailing temporal pulse in the 
SH (this often occurs close to or at phase matching in 
presence of GVM, see also [23[ ). This behaviour is also 
reflected in the SH spectrum, see (d) and cut in (f), which 
shows a very strong and extremely narrow red-shifted 
component building up, which eventually becomes the 
dominating contribution. As we discuss below its spec- 
tral position can accurately be predicted by the nonlocal 
theory that was recently developed by us 0, [H[ . We 
believe that this strong and long SH trailing component 
actually causes the trailing part of the FW to be strongly 
depleted, and that this is the main reason for the asym- 
metrical FW shape. 

The question is now: can we increase the effective soli- 
ton order and achieve further compression below 100 fs 
as to approach the limit predicted by Eq. ©? This turns 
out to be impossible: when N e g is increased the GVM 
effects become stronger, making the compressed pulse 
more distorted. This is clearly observed in Fig. [6l where 
we increase N e g and compare with the compression of 
Fig. O For N e g = 2.0 the compressed FW pulse in (a) is 
still quite short, but clearly is less clean. For N e g = 2.5 
the compressed pulse instead becomes quite distorted. 
It is also evident in the SH time plots that the trailing 
DC-like component increases with N e g, while the central 
part in all cases is a sub- 100 fs FWHM pulse. It is quite 
weak because most of the converted SH energy is fed into 
the DC-like part of the pulse, which is connected to the 
strong spectral peak in the SH spectrum. This spectral 



peak becomes stronger with increased N e g (not shown), 
but does not change position as it does not dependent on 

iVeff- 

In order to understand the spectral content of the 
different temporal components, the cross-correlation 
frequency-resolved optical gating (XFROG) method is 
useful. The spectral strength is given by [24] 



Sj(z,T,Q) 



dte int £j(z,t)£ gatc (t-T) 



(12) 



where £ ga te(i) is a properly chosen gating pulse. The 
spectrograms of the compressed pulses in Fig. [B] for 
Neg = 1.4 are shown in Fig. [7] The FW compressed 
pulse is slightly blue-shifted (around 2 THz), and the 
compressed part (located at T ~ 200 fs) shows a signifi- 
cantly broader spectrum. 

The SH spectrum is very particular: the part of the 
pulse that propagates with the FW group velocity (the 
"locked" part) shows a quite clean short pulse. This group 



velocity locking of the SH has been observed before 19 , 
23| and can be understood from the nonlocal theory 16|, 
18j : the SH has a component that is basically slaved to 
the FW due to the cascading nonlinearities. In frequency 
domain it can be compactly expressed as (l8l | 



U 2 (z,n) cx R-(n)F[Uf(z,t)} 



(13) 



where T[] denotes the forward Fourier transform, and Uj 
are properly normalized fields. Thus, the spectral content 
of the SH is slaved to the spectral content of the spectrum 
of Uf. The weight is provided by the nonlocal Raman- like 
response function in the nonstationary regime [l8j 



(n - Q-){Q - n+) 



(14) 



where f2± = f2 Q ±f2f,. These frequencies can be calculated 
(to 2. order) from the dispersion of the system as £l a = 

d 12 /k^ = -1.044 PHz and fl b = \2Ak/k 2 2} - Q.^ 1 / 2 = 
0.963 PHz. In the center around = 0, where F[U?(z, t)] 
is residing in this case, the response is quite flat: thus we 
get a SH component locked to the FW and when the FW 
compresses so does this SH component. 

Another striking feature of the SH spectrogram is the 
DC-like component: it is very evident as a long pulse cen- 
tered around £1 ~ —80 THz. Also this peak can be under- 
stood from Eq. (fl"3|) , because according to Eq. (fl4|) the 
nonlocal response function in the nonstationary regime 
has sharp resonance peaks in the response at = Q±. 
Inserting the dispersion values of the simulation we get 
S7+ = —81.6 THz in excellent correspondence with the 
observed peak position as the red dashed line indicates. 
Instead f2_ is located too far into the red side of the 
spectrum to affect the behaviour. 

Considering this spectral composition, it might even be 
possible to filter away the disturbing SH component at 
fl = f2 + , which in time-domain would give a quite decent 
SH pulse. In (c) we show that this is feasible: we pass the 
SH pulse through a super-Gaussian (n — 3) bandpass 
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FIG. 7. (Color online) XFROG-like spectrograms of the simulation in Fig. [5] at the optimal compression point z opt = 91 mm. 
The sech-shaped gating pulse had J|^ WHM = 70 fs, and the spectrograms are normalized to the peak value of Si. The top and 
side plots show the purely temporal and spectral traces, respectively, and are thus identical to Fig.[Sfe) and (f). The red dashed 
line in (b) indicates the value fi+ as calculated by the nonlocal theory. The spectrogram in (c) shows the SH passed through a 
3. order super-Gaussian bandpass filter centered at the SH carrier frequency A2 = 0.515 /im and with a FWHM of 100 THz. 



filter centered at u>2 and with a bandwidth of 100 THz 
FWHM (corresponding to 15 nm): this filters away the 
disturbing sharp peak, and a 80 fs FWHM pulse remains 
at A = 0.515 nm. The peak intensity in this short pulse 
is around 0.006/i n = 0.0414 GW/cm 2 . If we assume that 
it is created with 15 fi J pulse energy focused to wq = 0.5 
mm to achieve N e g = 1.4, and that the generated SH has 
roughly the same spot size, then the pulse energy of the 
filtered 80 fs pulse would be around 50 nJ. 

B. Case (2): 450 fs 50 fjj pulses 

In case (2) the pulse duration is longer, 450 fs. When 
the pulse duration is longer the soliton will for a fixed 
soliton order compress after a longer distance. This is 
because according to Eq. (JOj) z opt oc Ld,i oc T? n . How- 
ever, we may compensate for this by increasing the effec- 
tive soliton order enough to reach the limit governed by 
Eq. ([8]). For a 450 fs 50 /iJ pulse it is achieved around 
wo = 0.8 mm, see Fig. H(e), resulting in N c g ~ 2.0 — 2.5. 
This higher soliton order should make it possible to com- 




FIG. 8. (Color online) Numerical simulations using 450 fs 50 
/iJ input pulses and taking Ak — 45 mm^ 1 . The best pulse 
was observed for N c b = 2.0 (red curve) where pulse compres- 
sion occurs after 15 cm. The black curves show what happens 
as the effective soliton order increases (in which case the op- 
timal compression point occurs sooner). 



press in crystal lengths of around 10-15 cm, see Fig.UJc). 

In Fig. [5] we show some numerical simulations using 
these longer more energetic pulses. The best pulse ob- 
served shows a three-fold compression to Ai opt = 121 fs 
(FWHM) at iVeff = 2.0. The compression occurred after 
around 15 cm propagation, so spatial walk-off would be 
an issue here. Increasing the soliton order to N e s = 2.6 
the pulse becomes more distorted, but still compresses to 
around 150 fs FWHM after 9.5 cm, a more realistic in- 
teraction length. Finally, at N c g — 3.0 the pulse becomes 
too distorted as the GVM effects become stronger. 

In the two cases the pulses therefore eventually com- 
press to the same duration, which is the limit imposed 
by the nonlocal GVM effects. The more energetic pulses 
in case (2) allow for a more defocused pump beam so the 
compression should be less affected by diffraction. On the 
other hand, as the pulses are longer they compress later, 
so spatial walk-off is a more severe issue. A more optimal 
situation in both cases would therefore be more energetic 
pulses so the pump can be defocused with a factor 2-3. 
This would diminish spatial walk-off effects. 

IV. CONCLUSION 

Here we have shown that lithium niobate (LN) crystals 
in a type I cascaded SHG interaction can provide moder- 
ate compression of fs pulses from Yb-based fiber amplifier 
systems (1.03 /im wavelength). The phase mismatch was 
controlled through angle tuning (critical phase matching 
interaction). Using numerical simulations we found that 
the best compression was to around 120 fs FWHM after 
around 10 cm propagation. 

Better compression was prevented in part by strong 
GVM effects, caused by strong dispersion in the LN 
crystal, and competing material Kerr nonlinear effects. 
These are focusing of nature and counteract the defocus- 
ing Kerr-like nonlinearities from the cascaded SHG. In 
order to make the total nonlinear phase shift negative 
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the phase mismatch had to be taken quite low, and in 
this regime GVM effects dominate (the "nonstationary" 
regime). GVM imposes a strongly nonlocal temporal re- 
sponse in the cascaded nonlinearity that feeds most of 
the converted energy into a narrow red-shifted peak. In 
the temporal trace this gave a SH with a multi-ps long 
trailing component. The FW therefore experienced a dis- 
torted compression less the soliton order was kept very 
low. For such low soliton orders the compression distance 
increases substantially, but here the strong dispersion of 
the LN crystal actually becomes an advantage: due to a 
large GVD the soliton dynamics occur in much shorter 
crystals than usual, and the numerics indicated compres- 
sion in realistic crystal lengths (10 cm). 

It was noted that using low soliton orders gave a com- 
pressed pulse retaining most of the input pulse energy (in 
the cases we showed around 80%), and that the unavoid- 
able soliton pedistal was less pronounced. 

We also discussed the implications of using long crys- 
tals. Spatial walk-off will be an issue since it is a criti- 
cal phase matching scheme is used that exploits birefrin- 
gence, and also diffraction can be a problem. In order 
to counteract these detrimental effects the pump pulses 
need to be as energetic and short as possible. Two cases 
were highlighted taken from commercially available sys- 
tems, and we argued that diffraction should not prevent 
observing the predicted compression, but that some sort 
of walk-off compensation might be needed. Future sys- 
tems with more energetic pulses and reasonably short 
pulse durations (< 500 fs) would be able to beat the walk- 
off problem. Walk-off could also be prevented by using a 
noncritical type I phase matching scheme (9 = n/2) and 
increasing the temperature slightly to around 45° C. 

We finally noted that the peculiar SH shape in the non- 
stationary regime gave a very characteristic spectrogram: 
as mentioned above nonlocal GVM effects resulted in a 
sharp spectral red-shifted peak with a long multi-ps trail- 
ing temporal component. Another pulse component was 
instead locked to the group velocity of the compressed 
FW soliton. This locked visible pulse was located at the 
SH wavelength (515 nm), quite far from the red-shifted 
peak. We showed that a simple bandpass filter could ac- 
tually remove the detrimental red-shifted peak leaving 
a very clean 80 fs visible pulse (A = 515 nm). This is 
the opp osite approach compared to other studies, see e.g. 
[HHHl, where focus was on exploiting "spectral compres- 
sion" of fs pulses to obtain longer ps pulses. Despite that 
the cascaded SHG by nature has a low conversion effi- 
ciency, the pulse energy of this short visible pulse can 
easily be 50-100 nJ. Such pulses could be used for two- 
color ultra-fast energetic pump-probe spectroscopy. 

This study showed that cascaded quadratic pulse com- 
pression is possible even in a very dispersive nonlinear 
crystal. However, if compression occurs in a medium with 
stronger quadratic nonlinearities then it would be possi- 
ble to increase the phase mismatch, and thereby enter 
the stationary regime where the nonlocal GVM effects 
are much weaker. The benefit would be triple: cleaner 



compressed pulses could be generated, higher soliton or- 
ders could be used to achieve stronger compression, and 
it would occur in a shorter crystal. This conclusion is in 
line with what was noted previously in a fiber context 
[26| , where one of us found that the very dispersive na- 
ture of wave-guided cascaded SHG could be overcome if 
a strong enough quadratic nonlinearity is present. We are 
currently investigating other possible nonlinear crystals 
and phase matching conditions to achieve this. 
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Appendix A: LN crystal parameters 

LN is a negative uniaxial crystal of symmetry class 3m. 
Its low damage threshold due to photorefractive effects 
and problems with green induced IR absorption can be 
improved dramatically by doping the crystal, in partic- 
ular with MgO doping [23,[28j]. 1% MgO doping in sto- 
ichiometric LN (1% MgO:sLN) is enough to practically 
remove photorefractive effects and increase dramatically 
the damage threshold, while 5% is needed in congruent 
LN (5% MgOxLN) to do the same [28]. 1% MgO:sLN 
also has a shorter UV absorption edge (A = 0.31 /im). 

We here use 1% MgO:sLN, and the Selhneier equa- 
tions from [22|: note that for 1% MgO:sLN they only 
measured n e , but we checked that the 5% MgOxLN n 
Scllmcier equation matches well (at room temperature) 
the 1% sLN n a equation from (29[. The quadratic nonlin- 
ear coefficients have been measured at A = 1.06_um and 
are d 3 i — —4.7 pm/V and d 33 = 23.8 pm/V [3(|, while 
rf 2 2 = 2.1 pm/V [31| was measured for undoped LN. The 
fact that d,3id,22 < has been established in, e.g., [32[. 
The effective quadratic nonlinearity of the type I oo — > e 
interaction is given by Eq. Because d3ic?22 < [HJ 
the maximum nonlinearity is realized with <p = —tt/2. 

Appendix B: Anisotropic Kerr nonlinear refraction 

We previously studied type I cascaded SHG in a BBO 
crystal [16I - 4X8I ] . assuming an isotropic Kerr nonlinearity 

XeS,U - X e ff,22 ~ 6 X e S,12 

(3) 

where X e s jj are * ne ^" an< ^ ^ SPM coefficients, and 
12 ' ls the XPM coefficient. However, all quadratic non- 
linear crystals are anisotropic, and below we address this. 

Note first that the error made in assuming an isotropic 
response for the CQSC is probably small as the crucial 
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On the reduced form the cubic tensor becomes 



FIG. 9. Definition (in accordance with the IRE/IEEE stan- 
dard of the crystal coordinate system xyz relative to the 
beam propagation direction indicated by k. 



parameter is the FW SPM coefficient. As we will see now 
for type I this is identical in the isotropic and in the 
anisotropic cases. However, it should be emphasized that 
the various experimental attempts to measure the Kerr 
nonlinear refractive index of nonlinear crystals do not al- 
ways measure the tensor component relevant to our pur- 
pose, namely the en component, see Table U later. The 
analysis presented here should help understanding what 
exactly has been measured, and put the results into the 
context of cascaded quadratic soliton compression. 

For a nonlinear crystal in the symmetry group 3m (LN 
and BBO) there are 37 nonzero elements for the x'' 3 '' 
tensor, and of these only 14 are independent (34| 

xxxx = yyyy = xxyy + xyxy + xyyx 

xxzz = xzxz — xzzx = yyzz — yzyz = yzzy 

= zyyz = zyzy = zzyy = zxxz = zxzx — zzxx 
xxyy — xyxy = xyyx — yxxy — yxyx — yyxx 
xxyz = xxzy = xyxz = xyzx = xzxy = xzyx 

= ~yyy z = ~yyzy = ~y z yy = yxxz = yxzx 

= yzxx = —zyyy = zxxy — zxyx — zyxx 
zzzz (B2) 

where Kleinman symmetry has been invoked, and the po- 
larization relative to the crystal coordinate system is de- 
fined in Fig. [9l Under Kleinman symmetry the nonlinear 
coefficients are assumed dispersionless and the criterion 
for this assumption is that the system is far from any 

(3) 

resonances. Using the notation xljki = C A"« where 

for /i : x — * 1 y — > 2 z — > 3 

for m : xxx — > 1 yyy — > 2 zzz — > 3 yzz — » 4 
yyz — > 5 xzz — » 6 xxz — > 7 xyy — » 8 

xxy — > 9 xyz — > (B3) 

these tensor components are equivalent to 

en = c 22 = 3ci 8 

Ci6 = c 24 = c 35 = c 37 

Cl8 = C 2 9 

ClO = -C25 = C27 = -C32 = C39 

c 33 (B4) 



(B5) 



en 

CH 

-cio C33 







ci6 f cio 






Ci 6 -Cio Cio 

Cl 6 Cio 



£11 n 
3 u 



Cl6 



These results conform with the IRE /IEEE standard [35| • 
We now want to evaluate the cubic nonlinear response 
for a type I interaction. Using the notation from [l7| the 
cubic nonlinear polarization response is 



P^ = eoX (3) iEEE 



(B6) 



We have here only considered an instantaneous (elec- 
tronic) cubic nonlinear response [36| . Let us consider the 
type I SHG interaction where two ordinarily polarized 
FW photons are converted to an extraordinarily polar- 
ized SH photon (00 — * e). In the coordinate system ac- 
cording to the IRE/IEEE standard [H, see Fig. EJ the 
unit vectors for o-polarized and e-polarized light are 



(B7) 



where walk-off has been neglected. 

We then introduce slowly varying envelopes polarized 
along arbitrary directions 





— sin <j) 




— cos 9 


cos</> 
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COS(f) 
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— cos 9 


sin <f> 
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E(t) = Re[ui£"i(t)e 



-%UJ\t 



u 2 £ 2 (t)e 



(B8) 



where Uj is the unit polarization vector. For type I SHG 
we have Ui = e° and u 2 = e e . The nonlinear slowly 
varying polarization response 



NLV 

then becomes 

?(3) _ 



P- 



NL,i 



v (3) 



Si, (B9) 



where i,j = 1,2 and j ^ i. We have here only included 

phase-matched components and frequency-mixing terms 

where 2wi_— 0J2 = 0. The numerical prefactor | is the 

if-factor [33] for a third order nonlinear effect creating 

an intensity dependent refractive index with degenerate 

(3) 

frequencies, and the factor 2 on the XPM terms X c s ij 
stems from the fact that the if-factor for cross-phase 
modulation with non-degenerate frequencies is |. 

For calculating the cubic nonlinear coefficients, it is 
convenient to use an effective cubic nonlinearity [381 ] 



xfs = u d • x (3) :u a u 6 u c 



u d ■ c ■ u 



(3) 



(BIO) 



a, 6, c, d = 1,2. Here Ud is the unit vector of the field un- 
der consideration; thus, if we are interested in calculating 
the cubic nonlinear polarization for the FW [taking i = 1 
in Eq. (|B9[) ]. then = uj.. The other three unit vec- 
tors u Qj b jC are the unit vectors of each field appearing in 
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Eq. (|B6[) , and can in the case we are considering here be 
either u x or u 2 according to the identity (|B8|) . Most com- 
binations are not phase matched or have 2u>i — u>2 =/= 0, 
and are therefore not included in Eq. (|B9|) . The rank 4 
tensor on reduced form, as given by Eq. (|B5|) for LN, 



can be used to find the tensor product x( 3 ):u a Uf,u c as a 
simple matrix-vector product c • where 



,(3) 



Ly 

Ly 

L.i: 

Lj 



L 

L z 
L z 
U, 



yyy 



j yzy -r ^ Z yy 
~t~ -^zxz ~t~ -^zzx 

1 " L xzx -\- Lzxx 
L, 



J yyx 



J xxz 
L X yy 

Lxxy ~T~ J^xyx ~T~ uyxx 
Lxzy T L zx y ~\~ Lyxz ~\~ Ly 



Ln, 



+ L Z 



(Bll) 



Here Lju = u a ^Ub,kUc,i where the jkl indices refer to the 
x, y or z components of the unit vectors. 

It is convenient at this stage to simplify the notation 
based on the type I SHG interaction we are interested 
in. The effective cubic nonlinearity (|B10[) then reduces 
to the nonlinear coefficients appearing in Eq. (IB9[) 



(3) 
Acff ,ij 



•X (3) :u iU ,- 



(B12) 



The SPM terms can now be calculated as follows. The 
FW SPM interaction has i = j = 1 in Eq. (IB12|) . and is 
an ooo — > o process: Ui = e°. The SH SPM interaction 
has i — j = 2 and is an eee — > e process, so u 2 = e e . 

We then need to calculate x ■" a 4 u i u i using the reduced 
notation. Since for the SPM terms all the unit vectors in 
u' 3 ) are degenerate in frequency, all Lm components in a 



given vector entry are identical, e.g. L yzz 
We then get for the FW 



L, 



X (3) :uiuiui = 



— en sin ( 
Cn cos 4> 
-Cio cos 3 



(B13) 



A similar expression can be calculated for the SH SPM 
component, although it is substantially more complex. 
In the final step we carry out the vector dot product of 
these vectors with u i: as dictated by Eq. (IB12|) . and get 
for the FW (i = 1) and the SH (i = 2) 



(3) 

Xcff ,11 = Cll 



Y (3) _ 
Acff, 22 " 



-4cio sin 9 cos 3 9 sin 3^ + cu cos 
+|ci 6 sin 2 26> + c 33 sin 4 ^ 



(B14) 



(B15) 



For the XPM terms note that the three unit vec- 
tors used to calculate Eq. (|B1 1|) are non-degenerate in 

frequency. As an example, for x^^UiUi terms like 




FIG. 10. The calculated effective Kerr nonlinear contributions 
from cascaded SHG to the measured Kerr nonlinear refractive 
index by Kulagin et al. [3S| . The beam propagates with 9 = 
tt/2 into a LN crystal, and ip denotes the polarization angle 
{^p = gives o-polarized light, while ij) — tt/2 gives e-polarized 
light). The angle (j) wa s not reported, but we checked it has 
little influence on the ng HG value shown here. 
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L yyx must be evaluated, whose compo- 
= — cos 9 cos 3 (f> and L yxy = L yyx = 

This gives x& 12 = Xe$ 21 and S H3 



y (3) 

A c ff, 12 — 3 



1 en cos 2 9 + cie sin 2 9 + c w sin 29 sin 3^B16) 



The next step is to obtain the the values for LN of 
each component in Eqs. (|B14[) - (|B16|) . The value of the 
cubic nonlinear refractive index has been measured by 
many authors and for many different pulse durations and 
crystal cuts. In Tab.|T]the x*- 3 -* tensor components and the 
n Kcrr are reported in electrostatic units values, and the 
latter is also given in SI units (see App. IClfor details). 

In one of the earliest studies the tensorial nature of LN 
was studied Another early study found that c\\ = 
3cio [ill . Later studies used Z-scan methods and often a 
nonlinear refractive index value was found without any 
mentioning of the tensorial nature of the cubic nonlinear 
susceptibility. The cascaded quadratic contributions were 
also often forgotten or neglected. 

A recent study by Kulagin et al. went into a detailed 
experimental determination of the various cubic tensor 
components of LN, and found cu = 2.4 x 1CP 13 esu 
at A = 1.06 //m, and that cis = 1.2ci6 = ITC33 [39| . 
Through the relation en = 3cis, see Eq. (|B4|). the other 
coefficients are Ci6 = cn/3.6, C33 = cn/4.2. A problem 
with this study is that the cascaded quadratic nonlin- 
ear contributions to the observed Z-scan results were ne- 
glected. Instead, based on an analysis of the anisotropic 
Kerr tensor components the Z-scan transmission function 
was calculated, and the various tensor components were 
found by fitting to experimental data. In the experiment 
the pump propagated with 9 = tt/2, i.e. with the k- vector 
perpendicular to the OA. The angle of the polarization 
vector was then varied; this gives either pure o-polarized 
light, pure e-polarized light, or a linear mixture. 

We have done an analysis of the various cascaded SHG 
processes that come into play (00 — > o, 00 — > e, oe — > e, 
oe —> o, ee — > e, and ee — * o), evaluated their respective 
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A 

M 


(31 
Xcff 

[10~ 13 esu] 


J 

^Kcrr 

[10~ 13 esu] 


j 

n Kcn 

[10~ 20 m 2 /W] 


tFWHM 
[PS] 


Rep. 


e 

[deg] 


pol 


n 


Cij 


Ref. 


Note 


1064 


1.1 


48 


9.1 


30 


single 


90 


e 


2.2 


C33 


[40] 


:r-cut 


1064 


0.73 




6.0 


55 


2 Hz 


90 





2.2337^ 


Cll 


[41] 


paraxial fit 


1064 


0.66 




5.4 


55 


2 Hz 


90 





2.2337~ 


Cll 


[41] 


Gaussian fit 


1064 


2.4 


10 


19 


55 


2 Hz 


90 





2.2337 


Cll 


[39] 


Fit to transmission curve 


1064 


0.80 


3.4 


6.3 


55 


2 Hz 


90 


e+o 


2.2337 


Cl2,Clg 


[39] 


", C12 = cn/3 


1064 


0.57 


2.8 


4.9 


55 


2 Hz 


90 


e 


2.1495 


C33 


[39] 


", c 33 = C12/I.4 


1064 


0.67 


2.9 


5.5 


55 


2 Hz 


90 


e+o 


2.1912 


C23, Cl6 


[39] 


", C 23 = Cia/1.2 


800 


1.8 


7.8 


15 


0.42 


1 kHz 


? 


7 


2.1677 a 


7 


[42] 


£-CUt, Z-CUt 


lOU 


z.b 


11.0 


2U 


O.lo 


ID MHz 


U 





z.Zool 


Cll 


143] 


070 Mg(J:-LlN, z-cut 


577 


L6 


6.6 


12 


5,000 


40 Hz 








2.301 a 


ci8 


[44] 


C18 = Cn/3 


532 


10 


44 


83 


22 


single 


90 


e 


2.23 


C33 


[40] 


x-cut 


532 


6.6 


28 


53 


25 


10 Hz 








2.2244 a 


Cll 




z-cut 


520 


5.0 


21 


39 


0.2 


1 kHz 


90 


e 


2.24 


C33 


w] 


5% MgO 0.06% Fe cLN 



a Linear refractive index not provided; this value was calculated by us for conversion purposes. 



TABLE I. Nonlinear Kerr refractive index of LN measured mainly by the Z-scan method [13] • The underlined results are the 
values reported. The other entries have been calculated using Eqs. (1C8[) - (|C11[) . 



deff-values and phase mismatch values as the input po- 
larization angle changes. In total we arrived at a strongly 
varying cascaded contribution shown in Fig. 1101 At %jj = 
the contribution from ng HG is focusing, implying that the 
en component in Kulagin et al. might be too high with a 
factor of 7.0 x 10 -20 m 2 /W. There are also strongly defo- 
cusing contributions at other polarization angles, which 
should give rise to an underestimated value of the other 
tensor components. Moreover, the overall shape reminds 
strongly of the shape found in Fig. 5 in [39(: the focus- 
ing peaks from cascaded SHG could explain the valleys 
found there, and the defocusing valleys from cascaded 
SHG could instead explain the peaks. In summary we 
believe the Cn value to be too high, and the relation to 
the other tensor components to be dubious. 

There are other issues with the Z-scan method: If the 
repetition rate is too high, there will also be contributions 
to the measured n^ cl . r from thermal effects as well as two- 
photon excited free carriers [4£| , and hence n^eir does not 
contain just the instantaneous electronic response, as it is 
supposed to. Similarly conclusions can be made for pulses 
longer than 1 ps. For more on these issues, see e.g. [50j |. 

For the CQSC system the by far most important com- 
ponent is the FW SPM coefficient n^ err)11 . The SH SPM 
and the XPM coefficients only play minor roles in ex- 
treme cases close to transitions (e.g., close to the soliton 
existence line in Fig. [3]) . We checked in the cases we stud- 
ied in this paper that even increasing the SH SPM and 
XPM Kerr coefficient several times the isotropic values 
did not significantly change the compression results. 

Therefore until detailed reliable measurements of the 
cubic tensorial components of LN become available, we 
decided to use an isotropic Kerr response, and focus on 
using a realistic value of the FW SPM coefficient. The 
best choice seems to be nL = 20 x 10~ 20 m 2 /W at 

A = 0.78 /jm found in Ref. [431 ]. In this experiment they 

(3) 

have — and thus what they measure is x c s = c n- 
For orthogonal input polarization (corresponding to <p — 



0,7r/2, both cases o-polarized) they find the same value 

(3) 

as they should since this x c g does not depend on 0, cf. 
Eq. (|B14[) . Since they used fs pulses problems with long 
pulses are avoided. The high repetition rate could cause 
concern, but they checked that lowering it to below 1 
MHz did not change the results. Finally, the contribu- 
tion from the cascaded nonlinearities should be low: we 
estimate |ng HG | < 10~ 21 m 2 /W. 

As discussed later in App.lDlwe use Miller's rule to con- 
vert the nonlinear coefficients to the Ai = 1.03 /im that 
we use in the simulations in Sec. IIIII This implies that 
in the numerics we use n^ el . r n = 18.0 x 10~ 20 m 2 /W, 

n Kcrri2 = 6 -° x 10_2 ° m 2 /W, and < orr22 = 18.3 x 
10 -2d m 2/ w 

Appendix C: Conversion relations 

Often the nonlinear susceptibility is reported in Gaus- 
sian cgs units (esu) instead of the SI mks units. The con- 
version between esu and SI is 

X$ = 47T X ^(10 4 / C ) 2 (CI) 

where c is the speed of light in SI units. The 47r comes 
from the Gaussian unit definition of the electric displace- 
ment D = E + 47rP, and the 10 4 /c comes from converting 
statvolt/cm to V/m. 

In most cases the nonlinear Kerr refractive index is 
used. It is usually defined as the intensity-dependent 
change An in the refractive index observed by the light 

n = n + = n + rtKcrr^ol 2 = n + n^ err ^C2) 

Here no represents the linear refractive index, £q and Iq 
the input electric field and intensity, respectively. In our 
case the total polarization (linear and cubic, in absence of 

quadratic nonlinearities) can be written as Pj = + 

^NLi = £ o( £ i + £NL,i)£j- Now writing the sum of the 
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linear and nonlinear relative permittivities as £; + £nl.; = 
(n, + |An t ) 2 ~ n 2 + riiArii (here we take An; -C rij) 
then we can write the change in refractive index due to 
the Kerr nonlinearity on the form 



An; ~ nKerr.iil^il + 2n K ci-T,ij \£j I (C3) 

When comparing with Eq. (|B9[) we get in SI units [5l[ 

3 

4n, 



/'K,,r,i J -(SI) = £-xgy(SI), i,i = l,2 (C4) 



Note that the numerical prefactor 3/4 is the iT- factor dis- 
cussed above. Adopting the intensity notation the change 
in refractive index is Am ~ 2(n^ eTlii Ii+2n^ exTij Ij), and 
since in SI units Ii = ieo n i c l£i| 2 ; we get 



;(SI) 



1 



TljEoC 

3 



(SI) 



AniUneoc 



xS«(si) 



(C5) 
(C6) 



With Gaussian cgs units we would instead get [5l| 



7 



n 



Kerr,ij 



(esu) 
(esu) 



— Y (3) 
„ . A off ,ij 



(esu) 



4tt 



(esu) 



12ir 2 

nitijC 



xiiL( esu ) 



(C7) 
(C8) 

(C9) 



We have here used that in Gaussian units the intensity is 
ij(esu) = (8ir)~ 1 nic\£i(esu)\ 2 . The iGfactor appears also 
in Eq. ([C7]) as jAtt = 3ir. Note that c is still in SI units 
in these expressions. 

The connection between the Gaussian and SI systems 
can best be done via Eq. (|C1|) and (|C6[) to give [13, HH 



(3) 
A c ff,ij 

"Kerr ,ij 



TljC 



(SI) 



(esu) = j^n^ err!ii (SI) 



(CIO) 
(Cll) 



where we have used that the SI system defines £$c 2 = 
l/u = 10 7 /4tt using c = 299 792 458 m/s exactly. 

Note that often the definition of the Kerr nonlinear 
refractive index is n = no + An = no + n.Kcrr|£ | 2 = 
rio + njj err J (in Ref. [l7| we used this notation), which 
introduces an additional factor of 2 between nxcrr and 
while the relation between n£-„, and y*- 3 ) is unaf- 



'Kcrr' 



fected. Thus, working with and n^ eTT is the safest 
because one never has to worry about this factor of 2; 
as an example Eq. (IC10|) is still valid, while with the 
alternative definition Eq. (|C11|) becomes riKcrr(esu) = 
(n o c/807r)n^ crr (SI) ^ . 



Appendix D: Wavelength scaling of the nonlinear 
susceptibility: Miller's delta 

In the results presented here we account for the wave- 
length dependence of the nonlinear coefficients by using 
Miller's rule, which states that the following coefficients 
(the Miller's delta) are frequency independent [52| 



5^ = 



(2) 



(1) (1) (1) 
A-fci A.jj A-kk 



i,j,k = x,y,z 



(Dl) 



and we remind that the linear susceptibility is 1 + XiP = 
nf. A similar relation holds for the cubic nonlinearity 



£< 3 ) = 



(3) 
A-ijkl 



(1) (1) ( 1 ) v ( 1 )' 
Xii Xjj Xkk Xll 



i,j,k,l = x,y,z (D2) 



We remark that Miller's delta is based on an anharmonic 
oscillator with a single resonant frequency and only gives 
a ballpark estimate of the value, and thus is not to be 
expected to have a large accuracy (see, e.g., [H, H3]). 
However, it has been shown to work decently for most 
nonlinear crystals [55| . 
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